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ABSTRACT 



Aims. We study a sample composed of 28 of the brightest stars in the Arches cluster. Our aim is to constrain their stellar and wind 
properties and to establish their nature and evolutionary status. 

Methods. We analyze K-band spectra obtained with the integral field spectrograph SINFONI on the VLT. Atmosphere models com- 
puted with the code CMFGEN are used to derive the effective temperatures, luminosities, stellar abundances, mass loss rates and wind 
terminal velocities. 

Results. We find that the stars in our sample are either H-rich WN7-9 stars (WN7-9h) or supergiants, two being classified as OIf + . 
All stars are 2^1 Myr old. There is marginal evidence for a younger age among the most massive stars. The WN7-9h stars reach 
luminosities as large as 2 x 10 6 L o , consistent with initial masses of ~ 120 M Q . They are still quite H-rich, but show both N enhance- 
ment and C depletion. They are thus identified as core H-burning objects showing products of the CNO equilibrium at their surface. 
Their progenitors are most likely supergiants of spectral types earlier than 04-6 and initial masses > 60 M Q . Their winds follow a 
well defined modified wind momentum - luminosity relation (WLR): this is a strong indication that they are radiatively driven. Stellar 
abundances tend to favor a slightly super solar metallicity, at least for the lightest metals. We note however that the evolutionary 
models seem to under-predict the degree of N enrichment. 
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1. Introduction 

The center of our Galaxy is a unique environment to study 
massive stars. It harbors three massive clusters - the Arches, 
Quintuplet and central cluster - which together contain about 
30% of the number of Wolf-Rayet stars known in the Galaxy 
(Ivan der H ucht 2006). Interestingly, the three clusters have dif- 
ferent ages, ranging from ~ 2 Myr for the Arches to ~ 6 Myr 
for the central cluster. Consequently, they host different popula- 
tions of massive stars and sample the entire upper part of the HR 
diagram. Studying their stellar content gives us a unique oppor- 
tunity to understand how massive stars evolve. 

Although there is a global framework explaining the evolu- 
tion of stars more massive than > 20 M H , a quantit a tive d escrip- 
tion is still lacking. According to Crowth er et al.l {1995), stars 
with masses in the range 25-60 M Q experience the sequence O 
-> Of -> LBV or RSG -> WN8 -> WNE -> WC, while for more 
massive stars, the sequence O — » Of — » WNL+ab s — > WN7 (— > 
WNE) — > WC is preferred. lLanger et all d 19941) favor another 
scenario in which all stars have a H-rich WN phase prior to a 
LBV event: O -> H-rich WN -> LBV -» H-poor WN -> H-free 
WN — > WC. We see that there are still some qualitative differ- 
ences between the proposed scenarios. Further, the evolutionary 
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sequences are not mature enough to allow a refinement of the 
classification of the stars in the different evolutionary states. For 
example, the spectral types of the O or WC stars entering the 
above scenarios are not specified. The question of whether or 
not all massive stars go through a LBV phase is also not an- 
swered. This is an important caveat, especially since this phase 
has recently been argued to be the one in w hich most of the mass 
removal happens (Smith & Owocki 2006|). 

In a previous study (M artins et alj l200"7t) . we analyzed 18 
massive stars in the central cluster of the Galaxy. This clus- 
ter is especiall y intriguing since i t hosts the superm assive black 
hole SgrA* dGenzel et al.lll996t iGhez et alj|1998l) . In spite of 
the drasti c tidal forces, several tens of massive stars formed 
recen tly dAllen et all 11990b iKrabbe et all 11995b iPaumard et ail 
2006). Some of them are ap proaching the black hole at distances 
of on ly a few light hours dGhez et al.l 120031: lEisenhaue r et al.l 
120051) . The presence of young stars in the Galactic Center 
together with the apparent implausibility of forming stars so 
close to the central supermassive black hole is a puzzle usu- 
ally referred to as "the parad ox of youth". Studyin g the dy- 
namics of these young stars. IPaumard et al] d2006l) (see also 
iLevin & Beloborodovl2003l:lGenzel et al.l2003l) have shown that 
they orbit SgrA* in two counter-rotating disks. Together with 
other evidences (total mass and structure of the disks), this points 
to a local, "in-situ" star formation event. The detailed analysis 
of the post-main sequence massive stars has revealed that, sur- 
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prisingly, their evolution follows almost perfe ctly the predictions 
of evolutionary models (Ma rtins et al.ll2007l) . This implies that 
whatever the exact formation mechanism is, the subsequent evo- 
lution is not different from that predicted for normal stars. We 
found that all stars seem to have progenitors in the mass range 
25-60 M and that t hey follow relat i vely w ell the evolutionary 
scenario proposed bv lCrowther et al.l {1995) for this mass range. 
We have been able to refine this scenario, pinpointing the rela- 
tion between different spectral types: O — > (Ofpe/WN9 <-> LBV) 

WN8 WN8/WC9 WC9. This was made possible by 
the detailed study of stellar abundances in various Wolf-Rayet 
stars and related objects. Abundance analysis is a powerful tool 
to constrain stellar evolution since it gives direct access to the 
evolutionary state of a star. 

The above study focused on stars in the mass range 25-60 
M due to the age of the central cluster (more massive stars do 
not exist any more). In order to constrain stellar evolution at very 
high mass, we need to study younger clusters. The Arches cluster 
in the Galactic Center, only 30 pc away from the central cluster, 
is the ideal target. Not only is it believed to be quite young (2 to 
4.5 Myr, see Figer et al. 1999, Blum et al. 2001 and the present 
study), but it also shares the same environment as the central 
cluster, and hence have the same metallicity. Hence, its study en- 
sures to obtain a homogeneous view of stellar evolution among 
all types of massive stars in the Galactic Ce nter. The Arches 
cluster , first discovered bv lNagata et al.l d 19951) and lCotera et alj 
(119961) . is also one the most massive and densest cluster of the 
Galaxy. iFiger et ail d 19991) first showed that the mass function 
(MF) of its central regions might be shallower th an the typ- 
ical S alpeter IMF. This result was confirmed by IStolte et aU 
(120021) . Although there are indication that hig h mass star forma- 
tion might be favored in th e Galactic Center dMorris & Ser abvn 
ll996tlKlessen et al ] |2007l) . recent simu l ations of the dynamical 
evolution of the cluster by iKim et ail d2006l) indicate that we 
might in fact witness the effects of mass segregation in the clus- 
ter core rather than an actual top-heavy initial mass function. 
Whatever the physical reason, there are nearly 100 massive stars 
in the Arches cluster. From K-b and s pectroscopy of the bright- 
est members. TBlum et alj d200ll) and Figer et alj d2002l) identi- 
fied several late WN sta rs and early O super giants. The analysis 
of five of these stars bv lNajarro et al.l (|2004) revealed, by an in- 
direct method, that their metallicity was close to solar. 

Here, we analyze a much larger sample (28 stars in total) in 
order to better constrain their stellar and wind parameters. We 
rely on high quality data obtained with the integral field spec- 
trograph SINFONI on the VLT. The K-band spectra extracted 
from this data set are analyze d with atmosphe r e mo dels com- 
puted with the code CMFGEN lffillier & Miller! ([19981) . In Sect. 
|2]we describe the observations, our sample and the spectroscopic 
classification; in Sect.[3]our method to analyze the stars are pre- 
sented; the results are summarized in Sect. [4] and discussed in 
Sects. |5]|6]and|7] We give our conclusions in Sect. [8] 

2. Spectroscopic data 

2.1. Observations and data reduction 

The Arches cluster was observed in service mode between 
May 3 rd and June 27 th 2005 with SINFON I on the ESO/VLT 
dEisenhauer et all 120031 iBonnet et ail 12004). K band data were 
obtained with a pixel scale of 0.1". Adaptive optics was used 
to improve the spatial resolution. The seeing varied between 
0.5 and 1.2 arcseconds during the different runs. Four sub-fields 
were observed at the core of the cluster, as well as 12 fields in the 



outer part. For each sub-field, the integration time on source was 
240 seconds. One sky exposure was obtained every two object 
exposures. Early B stars were observed as telluric standards. 

Data w ere reduced with t he SPR ED software dAbuter et al.l 

12006b as in lEisenhauer et all (120051) and lPaumard et alj d2006). 
The reduction steps include: sky subtraction, flat field and bad 
pixel correction, distortion correction, wavelength calibration 
and atmosphere correction. In the last step, the telluric standard 
is used after its intrinsic Bry line is removed by a simple interpo- 
lation of the continuum red and blueward of the line. Individual 
frames were subsequently combined to obtain mosaics of the 
observed reg i ons (w hen frames overlap). We refer the reader to 
lAbuter et alj (120061) for a comprehensive description of the soft- 
ware used. Spectra were then carefully extracted by selection 
of individual "source" pixels from which "background" pixels 
are removed to correct for light contamination. The final spectra 
have a resolution of ~ 4000 and a signal to noise ratio of 10 to 
80 depending on the brightness of the star. 

2.2. Sample 

We selected the stars with high enough signal to noise ratio 
spectra (S/N > 10) for a subsequent quantitative analysis with 
atmosphere models. Equivalently, this means that we studied 
the brightest members of the Arches cluster. The list is given 
in Table [0 The name of the stars is taken from the list of 
lFigeretalTd2002l). NICMOS ph otometry in the F205W filter 
was taken from lFiger et alJ d2002). and was assumed to be equiv- 
alent to K-band p hotometry. We also inclu ded star number 1 of 
iBlum et alj d200ll) which is not in the list o flFiger et ail d2002l) . It 
is designated by the name B 1 . Its K-band m agnitude is taken as 
the 2.14yum magnitude of IBlum et al.l d200ll) . This wavelength 
range is free of strong line. We estimated the absolute K-band 
mag nitudes adopting a dist ance to the Galactic Center of 7.62 
kpc (Eisenhauer et al. 2005). We also adopted a constant extinc- 
tion Ak = 2.8 for all stars. This value is slig htly l ower than 
the av erage A K derived by IStolte et ail d2002l) and IKim et all 
(120061) . However, as noticed by these two studies, the extinc- 
tion is smaller in the clu ster center ( i nner 5 "). In this region, Ak 
is between 2.6 and 2.95 Stol te et alj d2002l) . This behavior is in- 
terpreted as evidence for swiping of dust by stellar winds and/or 
photo-evaporation by the intense UV radiation of massive stars. 
Since most of the stars in our selected sample are in the cluster 
center, the choice of Ak = 2.8 is a reasonable assumption. From 
the adopted distance and extinction, we can derive the absolute 
magnitudes of our sample stars: they are reported in column 3 of 
TableQ] 

2.3. Spectral classification 

Spectral classification in the K band is more difficult than in 
the classical optical range due to the limited number of lines. 
However, catalogs of K band spectra of objects with spec- 
tral types derived from optical studies are becoming available , 
makin g the spectral classification easier ([M orris et al. 1 )19961: 
lHanson et ai1ll996t IFiger et al.1119971 lHanson et alj|2005l) . The 
main lines observed in our SINFONI spectra are the following: 
He i 2.058 [im, Hen 2.189 /urn, Hen 2.037 /urn, Hen 2.346 //m, 
Bry, Nin 2.247, 2.251 /urn, Civ 2.070-2.084 jum, Siiv 2.428 fim 
and the complex at 2.112-2.115 /mi (composed of Hei, Nm, 
C in and O m). They are identified in Fig. IA.ll to IA.4I The most 
prominent line, Bry, goes from a strong emission in the bright- 
est stars to an absorption profile when the stars become fainter. 
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Hen 2.189 fim shows a similar behavior. Hen 2.037 fim and 
Hen 2.346 fim have a weak P-Cygni profile when present. Civ 
2.070-2.084 fim, Nm 2.247, 2.251 fim and Si iv 2.428 fim, when 
present, are always in emission. Finally, the 2.112-2.115 fim 
complex is in emission in most of the spectra. 

The K band spectra we obtained are typical of late WN 
(WNL ) and early O type stars, in agreement with Fig er et all 
d2002l) . WN stars earlier than WN7 have He n 2. 189 fim stronger 
than Bry, which is observed in none of our sample stars. The 
distinction between late WN subclasses is difficult when only 
the K band is available dMorris etalJ|!996|) . In WN7 stars, 
the Hen 2.189 fan emission is strong (although not as much 
as Bry). In later type stars, Hen 2.189 fim is much weaker. 
In general, this morphology is associated with a strong Hei 
2.058 fan emission. However, in H-rich late type WN stars 
(the so-c alled WNh stars) Hei 2.058 fim can be seen in ab- 
sorption dCrowther & Smithlll996l;ICrowther & Bohannanl 19971; 
Boha nnan & Crowfherf 1 19991) . which can be qualitatively un- 
derstood by a lowe r H e content. Accordin g to the atlas of 
lHanson et all d 19961) and lHanson et all d2005l) . late type O stars 
have a weak Hen 2.189 fim absorption, and the 2.112-2.115 fim 
line complex in absorption or weak P-Cygni. Our sample stars 
do not contain these spectral morphologies. We thus have only 
early type O stars (04-6). Among them, main sequence stars can 
be distinguished from supergiants by the shape of Bry: it is in ab- 
sorption on the main sequence and either absent (because filled 
by wind emission) or in emission (usually with a weak central 
absorption) in supergiants. Some O supergiants have stronger 
lines than standard 04-61 st ars (especi a lly Br y) and are iden- 
tified as OIf + supergiantfl Figer et all d2002l) argued that the 
distinction between WNL (especially WN7) and OIf + is very 
difficult. Both types have strong emission lines, but WNL stars 
usually have Hen 2.189 fim in emission or at least with a P- 
Cygni profile while OIf + have ma inly He 11 2. 189 fim in absorp- 
tion. However. IConti et"aTI d 1995b showed that exceptions exist: 
HP 16691 and HP 190429 a re two OIf + with Hen 2.189 fim in 
emission. Fig er et ail d2002l) suggested that the presence of the 
Nm 2.247, 2.251 fim emission in WNL stars but not in OIf + 
supergiants could be used to break the degeneracy in the spec- 
troscopic classification. Based on these considerations, we clas- 
sify the Arches SINFONI spectra using the following criteria 
(restricted to early O and late WN stars): 

• WNL stars have a strong Bry and 2. 1 12-2. 115 fim emission. 
Bry is stronger than the 2.112-2.115 fim complex. They 
show N in 2.247, 2.25 1 fim (emission). He 1 2.058 fim is weak 
and/or in absorption, so all our WNL stars are H-rich (which 
is confirmed by our quantitative analysis, see below) and 
thus are classified as WNh. When Hen 2.189 fim is purely 
in emission the spectral type WN7-8h is assigned. When it 
is in absorption or shows a P-Cygni profile, the star is classi- 
fied as WN8-9h. 

• OIf + have Hen 2.189 fim in absorption, a weak or no Nm 
2.247, 2.251 fim line, the 2.112-2.115 fim line complex and 
Bry in emission. The strength of Bry is similar to that of the 
2.1 12-2.1 15 fim complex. 

• O supergiants (OI) have the same morphological properties 
as OIf + stars, except that Bry is in emission but weaker than 
the 2.1 12-2.1 15 fim complex, or in weak absorption. 



Table 1. Photometry of the stars analyzed in the pr esent paper. 
The st ars are identified by their number in the list of Fig er et aTl 
d2002l) . Observed magnitude s are also from this sou rce. A dis- 
tance of 7.62kpc is assumed ( Eisenhauer et al. 2 005b. as well a s 
a uniform extinction Ak = 2.8 in the K band (Stolte et al. 2002). 



Star 


ST 


m K 


M K 


Bl 


WN8-9h 


11.11 


-6.10 


Fl 


WN8-9h 


10.45 


-6.76 


F2 


WN8-9h 


11.18 


-6.03 


F3 


WN8-9h 


10.46 


-6.75 


F4 


WN7-8h 


10.37 


-6.84 


F5 


WN8-9h 


10.86 


-6.35 


F6 


WN8-9h 


10.37 


-6.84 


F7 


WN8-9h 


10.48 


-6.73 


F8 


WN8-9h 


10.76 


-6.45 


F9 


WN 8-9h 


10.77 


-6.44 


F10 


04-6If + 


11.46 


-5.75 


F12 


WN7-8h 


10.99 


-6.22 


F14 


WN8-9h 


11.22 


-5.99 


F15 


04-6If + 


11.27 


-5.94 


F16 


WN8-9h 


11.40 


-5.81 


F18 


04-61 


11.63 


-5.58 


F20 


04-61 


12.16 


-5.05 


F21 


04-61 


11.77 


-5.44 


F22 


04-61 


12.02 


-5.19 


F23 


04-61 


12.19 


-5.02 


F26 


04-61 


12.34 


-4.87 


F28 


04-61 


12.17 


-5.04 


F29 


04-61 


12.26 


-4.95 


F32 


04-61 


12.42 


-4.79 


F33 


04-61 


12.42 


-4.79 


F34 


04-61 


12.49 


-4.72 


F35 


04-61 


12.37 


-4.84 


F40 


04-61 


12.67 


-4.54 



1 Historically, the index f denotes stars with strong N in 4634-4641 A 
and He 11 4686 A emission, while the symbol "+" refers to stars showing 
Si iv 4089-41 16 A in emission. 



These criteria remain qualitative on purpose, so that broad 
groups of stars can be defined without preventing the possibil- 
ity that some stars are intermediate between the groups, as is 
likely to be the case in a population of massive stars with such a 
narrow age spread. The criteria we defined should also be seen 
as relative criteria to compare the different stars of our sample. 
Our final spectral classification for each star is given in Table [T] 
The spectra of our sample stars are also displayed in Fig. IA.1I 
[ASlOandlAll 



3. Spectroscopic modeling of individual stars 

In this section we describe our method to derive the stellar and 
wind properties of the selected stars. We also present an estimate 
of the uncertainties on the derived parameters. 



3.1. Atmosphere models 

We used the atmosphere code CMFGEN dHillier & Milled! 19981) 
to derive the stellar and wind properties of a sample of Wolf- 
Rayet and O stars. CMFGEN calculates non-LTE atmosphere 
models with winds and includes a robust treatment of line- 
blanketing. A d e tailed description of the code was given by 
Hillier & Miller (1998). Its main character istics are also pre- 
sented in Martins e tall d2004 l2005bl 12007b . Here, we simply 
highlight a few important features: 
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hydrodynamic structure: CMFGEN does not compute self- 
consistently the density (and velocity) structure of the at- 
mosphere. The standard procedure consists in adopting a 
pseudo-hydrostatic structure on top of which a "J3 velocity 
law" is connect ed. Such a law is exp ected from theoreti- 
cal ground (e.g. Pauldr ach et al. 1986b. Here, we used the 
TLUSTY structures of lLanz & Huberrvl d2003h for the hy- 
drostatic part of the atmosphe re. We chose § = 0.8 since it is 
the ty pical value for O stars dPuls et al.lll996tlRepolust et al.l 
2004). It also leads to good fits of the of the observed line 
profile. Only for stars Bl and F5 did we have to use val- 
ues of 1.2 and 1.8 respectively to better fit the overall shape 
of the emission lines, especially Bry: larger £1 produce more 
centrally peaked lines (for the adopted clumping factor, see 
below). High values o f B are also found for O supergiants 
dCrowther et al.ll2002at IffiHier et al.ll2003l) . The use of such 
a structure may be questionable in the case of extreme su- 
pergiants and Wolf-Rayet stars. However, the wind density 
in these stars is usually so large that the photosphere is be- 
yond the hydrostatic layers, so that the underlying structure 
has little impact on the observed spectrum. 
line-blanketing and super-levels: CMFGEN includes line- 
blanketing through the super-level approximation. In prac- 
tice, levels of similar energies are grouped in a single super- 
level which is then used to compute the atmospheric struc- 
ture. Within a super-level, individual levels have the same 
departure coefficient from LTE. This is a very convenient 
way to treat directly line-blanketing, without using statisti- 
cal methods such as opacity sampling. Usually, only levels 
with large energies are grouped into super-levels. This may 
affect the strength of infrared lines since they mainly arise 
from transitions between such high energy levels. As a con- 
sequence, we decided not to use super-levels for H, He, NIII 
and CIV which contribute most of the K band lines: we used 
the full atom. In addition to H, He, C and N, we included 
O, Si, S and Fe in our models. No other elements could be 
treated due to the increased memory requirement implied 
by the large number of levels (typically about 3000 levels 
and 1300 super-levels). The elements included are anyway 
the ones responsible for most of the line-blanketing effects. 
Tests with additional species confirmed that adding other el- 
ements did not significantly affect the results. Similarly, tests 
run with different super-levels assignments lead to minor 
changes insufficient to modify quantitatively our results. 
microturbulent velocity: a value of 15 km s _1 was used in the 
computation of the atmospheric structure. For the detailed 
emergent spectrum resulting from the formal solution of the 
radiative transfer equation, we adopt ed v tur b = 10 km s~ 



This i s a reasonable value for O stars (Villamariz & Herrerol 
2000). For WR stars, it might be a little too low. However, 
test models with v tur b = 50 km syndicate barely any change 
in the resulting line profiles, because they are dominated by 
the wind. 

clumping: CMFGEN allows the inclusion of clumping. In 
practice, a volume filling factor approach is used. The fill- 
ing factor is described by an exponential law starting from 
a value of 1.0 at the base of the wind and declining to 
/ in the outer atmosphere, where the velocity reaches v^. 
Given the limited number of diagnostics in the K band, 
we have adopted / = 0.1 in our computations. This is a 
standard value for WR stars dHamann & Koesterkd]1998 ; 
Hillier & Mme^ll999tlMorris et al.ll2000tlHilher et al.l2001 ; 
Crowther et al.l l2002al 120061) . For O stars, the amount of 
clumping is still a matter of debate. Values of / as low as 



0.01 have been derived b y lBouret et al.l d2005l) ; iMartins et ail 
d2005bl) : iFullerton et alT(l2006h. But l arger values are also 
fo und: / = 0.2 by iRepolust et al.l d2004l) or / = 0.1 
bvlCrowther et al.l d2002bl) : lMarfins et al l J 2005b). Recently, 
iPuls et all d2006l) also showed that the clumping factor might 
vary non monotonically with radius in the wind, contrary to 
the usual assumption of atmosphere models. Given these un- 
certainties, the adopted value for / (0.1) is not unrealistic. If 
the clumping factor was smaller, then our mass loss rates for 
O supergiants would be overestimated by JQ.l/f. 

3.2. Method 

We briefly describe here the method we used to constrain the 
main stellar and wind parameters. 

o Effective temperature: T e g was constrained from the strength 
of He i and He n lines, as in most studies of massive stars. In 
practice, we used Hei 2.058 pm, Hei 2.112 pm, Hen 2.189 
pm, Hen 2.037 pm and Hen 2.346 pm as the main diag- 
nostics. We note that T e g is defined at the radius in the at- 
mosphere model where the Rosseland optical depth reaches 
2/3. For comparison with evolutionary models, it is usually 
useful to define T* as the temperature where the opacity is 
20. This corresponds to a deeper, quasi hydrostatic layer of 
the atmosphere, which is more similar to the outer radius of 
the evolutionary models. In general, T e ff and T* are almost 
identical for O stars, while for WR stars with denser winds 
they can differ by several thousands of degrees. 

o Luminosity: the K-band flux was used as the main indicator. 
The luminosity was adjusted so that the K-band magnitude 
of the models could match the absolute magnitude of the 
stars. In practice, the K-band flux depends not only on the 
stellar luminosity, but also on the mass loss rate since in the 
near infrared the ionized wind produces free-free emission. 
Hence, the luminosity was derived in parallel to M. 

o Mass loss rate and He abundance: M and the ratio of H to He 
content were constrained from the strength of the emission 
of Bry and the He lines. A change of M leads to a general 
increase of the emission in all the lines, while an increase of 
He/H strengthens the He lines and weakens Bry. 

o C and N abundances: The carbon and nitrogen content was 
derived from Civ 2.070-2.084 pm and Nra 2.247, 2.251 pm 
respectively. We also use the Siiv 2.428 pm line when ob- 
served to constrain the Si abundance. We note however that 
this last abundance determination is less reliable than for the 
other elements since the spectrum is much noisier at the posi- 
tion of the Si iv 2.428 pm line (at the red end of the K band). 

o Terminal velocity: the terminal velocity of the wind (v M ) was 
determined from the width of the Bry and the extend of the 
absorption dip of the P-Cygni profile of He 1 2.058 pm (when 
present). When none of these indi cators could be use d, we 
simply adopted v„ = 2.6 x v esc dLamers et alJll995h . v esc 
being the escape velocity. 

Due to the absence of strong gravity indicators, we adopted 
logg = 3.25 for the coolest stars, and \ogg = 3.50 for the 
hottest. This is a reasonable a ssumption in view of the calibra- 
tions^ of Martin s et alj d2 005a). Finally, the solar abundances of 
iGreves se & Sauval ( 1998) were used as reference^ 



2 The solar abundances ha ve been recently revised for the lightest el- 
ements ( Asplund et al. 2004), but not for the heaviest. Hence, we prefer 
to stick to the old values. 
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3.3. Accuracy of parameters determination 

The determination of the stellar and wind parameters is a long it- 
erative process: most spectral diagnostics depend on several pa- 
rameters. To estimate the uncertainties on our determinations, 
one would ideally need to run tens of models covering the pa- 
rameter space from which some kind of chi-square procedure 
(to be defined) could provide statistical errors. In practice, this 
is not possible in the present approach since this would lead to a 
prohibitively long process: a model and the associated spectrum 
require between 24 and 48 hours of cpu time; sampling correctly 
the parameter space (10 to 20 models for each parameter, and 
~ 10 parameters) for each of the 28 stars would imply several 
months of computations. This is in addition to the time needed 
to actually find the best fit model for each star. 

Hence, we prefer to rely on a more empirical way to estimate 
the uncertainties. For this, we chose to run a few test models 
for two typical stars (one WN8-9h and one O supergiant). We 
varied the parameters around the values of the best fit model and 
judged by eye when the resulting spectrum was not satisfactory 
any more: this was used to define our uncertainty. This gives a 
reasonable estimate of the accuracy with which the stellar and 
wind parameters are derived. Fig.[T]illustrates this procedure for 
M and N/H. 

In practice, we focused on the models for star F2 (WN8-9h) 
and star F21 (04-61). They were chosen as typical of their class 
of objects, both in terms of spectral morphology and derived pa- 
rameters. The typical errors are: +3000K (2000K) on T e g for O 
stars (WNLh), +0.2 dex on log f- , 0.2 (0.1) logM, 100 km s" 1 
on terminal velocities and ±50% (±30%) on abundances. The 
uncertainty on log j- depends mainly on the adopted distance 
and extinction, and is thus similar for WNLh and O stars. The 
uncertainty on abundances does not take into account any pos- 
sible systematics due to uncertainties in atomic data and model 
assumptions. 

4. Derived stellar and wind parameters 

The results of our analysis are summarized in Table [2] The final 
fits are shown in Figs. IA.lllA.2llA3l and lA.4l A summary of the 
general properties is given in Sect. 14.11 while a comparison to 
previous studies is made in Sect. 14.21 

4.1. General properties 

To avoid lengthy discussions, we do not describe the stars one 
by one but focus on their global properties. They can be summa- 
rized as follows: 

o The stars have very similar T e g: out of 28 stars, 17 have 
32500 < r e ff < 37500. The average effective temperatures 
of WN7-9h and 04-61 stars are 33600 K and 37600K re- 
spectively. Considering the error on the temperature estimate 
(±2000-3000 K), there is very little dispersion in T e ff among 
our sample. This explains partly the similarity in the mor- 
phology of the K-band spectra. This has the important con- 
sequence that all stars lie almost on a vertical line in the HR 
diagram (see Sect.|5]l. Such a feature can be used to trace the 
age of the cluster. One can also note that the average tem- 
perature of the Q4-6 supergian ts is similar to the calibrated 
values of lMartins et alj d2005al) for 05-5.51 stars. 

o The stars are luminous: all the stars we analyzed have 
log 5.7, some of them reaching log j^- 6.35. WN7- 

9h stars have, on average, log^= 6.14, 04-6 supergiants 
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Fig. 1. Uncertainty estimate for M and N/H for one WNLh star 
(F2) and one O supergiant (F21). The solid line is the observed 
spectrum, while the colored broken lines are the models. The red 
dot-dashed lines correspond to the best fit model. In the other 
models, only M (left panels) or N/H (right panels) have been 
varied. 

have log^ = 5.84. Contrary to the case of effective tem- 
peratures, there is a clear luminosity difference between 
the classes of objects, WN7-9h stars being the most lu- 
minous. Luminosities in excess of 1O 6 L are indicative of 
very high mass stars. This characteristic will be further dis- 
cussed in Sect. [6] Early type supergiants/giants have lumi- 
n osities in agreement w ith the expectations: the calibrations 
of lMartins etail (l2005al) indicate log ^= 5.78-5.95 for 04- 
61 stars. 

o The stars are usually H rich: all the stars have a ratio He/H 
< 1.0. The WN7-9h stars are He enriched, but not exces- 
sively, with He/H between 0. 1 and 1 .0. Actually, only three 
WN7-9h stars have He/H > 0.5: all the other WNLh stars 
have a lower He content. It is interesting to note that some 
WNLh stars do not seem to be He enriched (stars Bl, Fl, F9, 
F14, F16). All the other types of stars show a solar helium 
abundance (with the exception of star F10 which is slightly 
enriched - He/H = 0.2). 

o The stars show various degrees of chemical enrichment: as 
expected, WNLh stars have the strongest nitrogen enrich- 
ment, with X(N) between 0.005 and 0.028. O supergiants, 
including the 04-6If + stars, have solar or slightly enriched 
N abundances (X(N) = 0.002-0.007). In parallel, WN stars 
are carbon deficient (X(C) < 1.3x 10 -3 ; exception: starF16) 
compared to O supergiants. This is a clear indication of CNO 
processing. 

o Mass loss rates are stronger in WN stars: all WN7-9h stars 
have M> \0 5 M Q yr" 1 (exception: star F16, M- 6.3 x 10~ 6 
M yr~'). The 04-6 supergiants have lower mass loss rates 
(in the range 2 - 4.5 x 10" 6 M Q yr -1 ). 

From this general overview, one sees that the stars analyzed 
here share some common properties (T e ff) but also are very dif- 



Table 2. Derived stellar and wind parameters. The typical errors are: ±3000 K on temperatures, ±0.2 dex on log ^- and logM, 100 km s 1 on terminal velocities and ±30% 
abundances. Terminal velocities are adopted from v ra = 2.6 x v esc for O stars (except the two 04-6If + supergiants). 
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ferent in terms of abundance patterns. This will be discussed in 
greater depth in Sect. 16.11 



4.2. Comparison to previous studies 

The only attempt to derive quantitative properties of the Arches 
massive stars is bv lNaiarro et al. (2004). The authors focused on 
five stars (F3, F4, F8, F10 and F15) and used a similar tech- 
nique to determine the stellar and wind parameters. They relied 
on Keck/NIRSPEC spectra with a high resolution (~ 23300) but 
a narrower spectral coverage than our VLT/SINFONI data (only 
four windows centered on He 1 2.058 pmi, He 1 2. 1 12 yrni, Bry and 
Nm 2.247, 2.251 fan were observed). 

The par ameters deriv e d in b oth the present study and the 
analysis of iNaiarro et alj d2004l) are summarized in Table [3] 
Generally, there is a rather good agreement between both stud- 
ies for the WN stars. Note in particular the similar luminosi- 
ties, mass loss rates, terminal velocities and He and N abun- 
dances. We find effective temperatures systematically larger (by 
~ 2000K), but the difference is within the uncertainties (except 
for star F4). The largest difference is found for the carbon abun- 
dance, although, with one exception (star F3), the discrepancy is 
only a factor of 2 or smaller. Here, we argue that INaiarro et al.l 
(2004) did not cover the full spectral range around 2.07/im to 
observe the Civ 2.070-2.084 ftm lines (see their Fig. 1), which 
is included in our SINFONI spectra. We have consequently a 
larger number of diagnostics and we are able to better derive the 
C content. 

The differences are larger for the two 04-6If + stars. We think 
the better quality of our spectra allows a better estimate of the 
terminal velocity (the full width of Bry is well observed), the C 
and N content (we unambiguously detect the C iv 2.070-2.084 
ymx lines and the Nm 2.247, 2.251 fim doublet). The mass loss 
rates being similar, the wind densities are not (due to the differ- 
ence in Voo), which partly explains the different luminosities. The 
different effective temperatures, well constrained by our well re- 
solved He i and Hen lines (especially Hen 2.189 ym) complete 
this explanation. 



5. HR diagram and cluster age 

The stars analyzed in the present study are placed in the HR 
diagram in Fig. [2] The Geneva evolut ionary tracks including ro- 
tation from Mevnet lfe Maederl d2005l) are used to build the dia- 
gram. Isochrones are also shown. As previously described, there 
is a clear difference in the position of stars of different spec- 
tral types. The WN7-9h stars (filled circles) are brighter than 
the normal O supergiants. The extreme supergiants are interme- 
diate. One can immediately conclude that the WN7-9h stars of 
the Arches cluster are very massive stars: only the 120 M evo- 
lutionary track is able to reach luminosities larger than 1O 6 L . 
Even the less luminous WN7-9h stars are accounted for only by 
the tracks with M > 6OM . One concludes that in the Arches 
cluster, the WN7-9h stars are the descendent of stars more mas- 
sive than 60 M . The position of the extreme early supergiants 
(the 04-6If + stars) overlaps with the position of the less lumi- 
nous WN7-9h stars. It is thus likely that they are closely related 
to them (see also the next sections). 

From the position of the stars in the HR diagram, one can at- 
tempt to estimate their age. For that, isochrones are indicated in 
Fig. [2] The most luminous WN7-9h stars are 2-3 Myr old. The 
O supergiants (with the exception of the 04-6If + stars) seem to 
span a slightly wider age range (2-4 Myr). There is however an 



overlap between the brightest supergiants and the faintest WN7- 
9h stars. This suggests that on average the most massive stars are 
slightly younger than less massive stars (in the mass range 40- 
120 M Q ): one can clearly exclude an age of 3 Myr for the most 
luminous WN7-9h stars, while some of the O supergiants with 
L ~ 1O 58 L could be ~ 4 Myr old. This may be an indication 
that the most massive stars formed at the end of the star forming 
event that gave birth to the Arches cluster. This would be con- 
sistent with the scenario according to which the most massive 
stars are the last to form in a starburst event since their presence 
immediately imply a strong negative feedback which removes 
material for star formation. We note however that we are prob- 
ing a small region of the HR diagram and this would need to be 
confirmed by a deeper study of intermediate and low mass stars. 
Within the uncertainty on the effective temperature and luminos- 
ity of the O supergiants, one cannot exclude either that they have 
the same age as the WNLh stars. Besides, binarity may change 
the picture. If the most luminous stars were found to be binary 
stars, the luminosity of each component would have to drop (by 
as much as 0.3 dex in case of equally luminous companions). 
This would translate into an older age for the stars, now closer 
to the 3Myr isochrones in Fig. [2] Hence, the suggestion of a late 
formation of the most massive stars still needs to be confirmed, 
but is worth being mentioned in view of the present results. We 
note also that in NGC3603, a galactic cluster quite similar to 
the Arches, there exists a population of pre- main sequence stars 
younger than the mos t massive components (IStolte et al.| [2006; 
lHaravama et al. 120071) . Clearly, a study of fainter components of 
the cluster is required to confirm the suggestion of a correlation 
bet ween age and initi al mass. 

Figer et ail (1 19991) derived an age of 2+1 Myr for the Arches 
cluster based on photome try of the massive components. A bet- 
ter estimate was given bv lBlum et ail (2001) who used informa- 
tion on the spectral types in combination to evolutionary models 
to constrain the age of the Arches to the range 2-4.5 Myr. Finally, 
comparing the types of Wolf-Rayet sta rs present in th e cluster to 
the predictions of starburst models of Mevnet ( 1995), as well as 
using the detailed properties of one star (F8). lFiger et al.1 (T2002) 
concluded that the age of the cluster is 2.5+0.5 Myr. Our deter- 
mination, which relies on a more quantitative basis and uses the 
most recent evolutionary tracks, agrees nicely with this estimate 
(especially if only WR stars are considered). 

Finally, one should mention that the above results were 
obtained using solar metallicity tracks. As we will see later, 
the metallicity in the cluster might be slightly super-solar. In 
that case, the age we derived would be only an upper limit. 
Indeed, for Z=2xZ , comparison to the corresponding evolu- 
tionary tracks show that all the stars we analyzed would be 
younger than 3 Myr. 



6. Nature of the most luminous stars 

In this section we investigate the nature of the Arches stars stud- 
ied here as well as the possible relation between the different 
types of stars (O, WNLh). We first focus on the chemical evolu- 
tionary status, then discuss the wind properties before drawing 
our conclusions. 



6.1. Chemical evolution 

In the previous section we have seen that the stars classified 
WN7-9h in our sample appear to be very massive stars. To bet- 
ter unravel the nature of these objects, the analysis of their abun- 
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Table 3. Comparison between our results and the study of Naiarro et al. (2004) for the five stars in common. For each star, the first 
row gives our results, and the second one the results of Naiarro et al. (2004). Only the main parameters are listed. 
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Fig. 2. HR diagram of the Arches cluster. Filled symbols are the 
stars analyzed in this work. The Geneva evoluti onary tracks in- 
cluding rotation from IMeyne t & Maeder (2005) are plotted, as 
well as isochrones. For clarity, only the first ~ 4 Myr of the tracks 
are plotted. The solid line in the right upper part of the figure is 
the Humphreys-Davidson limit. The typical error on the position 
of the stars is shown at the bottom left. 

dance pattern is a powerful tool since it informs about their evo- 
lutionary status. 

In Fig. [3] we plot the hydrogen mass fraction as a function 
of luminosity. The symbols have the same meaning as in Fig. [2] 
The relation from the Geneva evolutionary tracks are overplot- 
ted (solid lines). One can immediately conclude that the WN7- 
9h stars are the only ones of our sample showing H depletion. 
However, even some of the WN7-9h stars appear not to be He 
enriched (none of the O stars are H depleted). Overall, the H 
mass fraction of WN7-9h stars ranges between 0.2 and 0.7. We 
interpret this pattern as a sign that the Arches WNLh stars are 
objects which left the main sequence recently, some of them be- 
ing almost unevolved in terms of H depletion. This is a very im- 
portant conclusion, because it means that some (and maybe all) 




5.4 5.6 5.8 6 6.2 6.4 

log _L 



Fig. 3. Hydrogen mass fraction as a function of luminos- 
ity in the solar metalli city Geneva evolutionary models of 
IMeynet & Maeden d2005l) (solid line) and as derived in the stars 
analyzed in this work (symbols). 

of these stars are still core H burning objects. This implies that 
these stars are young, consistent with our age estimate (see Sect. 
0. The comparison with the evolutionary tracks in Fig. [3] also 
indicates that they have masses in the range 60-120 M Q , similar 
to what was inferred from their position in the HR diagram. 

A quantitative confirmation of the young evolutionary status 
of the Arches WNLh stars is given by Fig. [4] In this figure, the 
carbon mass fraction (X(C)) is shown as a function of the nitro- 
gen mass fraction (X(N)). According to stellar evolution, X(C) 
decreases in the earliest phases while X(N) increases when H 
is burnt through the CNO cycle. Then when He burning starts, 
C is produced at the expense of N (and He). Fig. |4]reveals that 
the WN7-9h stars are all carbon poor and N rich compared to 
the O stars of the sample (with the exception of star 14 which is 
more similar to O stars). In fact, they cluster in a rather small re- 
gion of the X(C)-X(N) diagram, which suggests that they show 
the pattern of CNO equilibrium. Contrarily to Fig. [3] there is a 
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0.004 




0.01 0.02 
X(N) 

Fig. 4. Carbon mass fraction as a function of Nitrogen mass frac- 
tion in the Geneva evolutionary models of Mevnet & Maeder 
(2005) (solid line: solar metallicity; dashed line: twice solar 
metallicity) and as derived in the stars analyzed in this work 
(symbols). The evolutionary tracks plotted are for 20, 40 and 
85 M Q . 



clear distinction between WN7-9h and O stars which are more C 
rich and N poor. This allows us to unambiguously state that the 
WN7-9h stars are evolved objects but are still at the beginning 
of their post main sequence evolution. 

Fig.|4]shows that although the trend of decreasing X(C) with 
increasing X(N) predicted by the evolutionary models is qualita- 
tively reproduced by the observations, quantitatively the agree- 
ment is not perfect. In particular, the WN7-9h stars span a wider 
range in X(N) (0.005 to 0.028) than expected from the models. In 
the models, at the level of X(C) seen in the WN7-9h stars, X(N) 
should be around 0.013 according to the solar metallicity track. 
The track at twice the solar metal content allows X(N) as large as 
0.025. But for the level of X(N) observed, the C content in this 
track should be up to 4-5 times larger. An explanation purely 
by a non solar metallicity is not satisfying. Different initial ro- 
tational velocities cannot fully explain this scatter in X(N) val- 
ues. Between two stars with initial rotational velocities of and 
300 km/s, the difference of maximum value of the N mass frac- 
tion X(N) is 0.002. This is about one tenth of the range spanned 
by the WNLh stars (from ~ 0.01 to -0.03, excluding star F14). 
Hence, a spread in initial rotational velocities cannot be fully 
responsible for the observed scatter in X(N). 

One can also speculate that the evolutionary models do not 
predict a strong enough N enrichment in the early phases of 
evolution of massive stars. One important ingredient which is 
still negle cted in most evolu t ionary models is magnetic fields. 
Recently, Ma eder & Mey net (2005) have shown that the pres- 
ence of magnetic fields in massive stars could favor solid body 
rotation and consequently chemical mixing. Their Fig. 10 re- 
veals that He and N abundances can be significantly increased 
compared to non magnetic models. Interestingly, the effect of 
magnetic fields on chemical enrichment are larger for older stars. 



This is to be compared to the larger spread in X(N) for more 
evolved stars (WNLh compared to O supergiants) in Fig. [4] We 
do not claim that magnetic field can explain all the trends seen 
in this figure, but it might be an important ingredient. In that 
context, it is worth noting that Trundle & LennorJ d2005l) report 
a similar discrepancy between derived (by atmosphere model- 
ing) and predicted (by evolutionary tracks) N abundances for B 
supergiants in the SMC. 

To further investigate the chemical evolutionary status of the 
Arches stars, we have plotted in Fig. [8] the ratio of N to He 
abundance (by mass) as a function of effective temperature. Here 
again, the Geneva evolutionary tracks for solar and twice solar 
metallicity are overplotted. Theoretically, the N/He ratio probes 
the first phases of evolution. Indeed, while both N and He are 
produced during the H burning phase, the relative increase of 
the abundances relative to the initial values is larger for N than 
for He, simply because He is already a main element in the star 
while N is not. In practice, the N mass fraction changes by an or- 
der of magnitude, while the He mass fractio n increases by only 
a fact or of 2-3 (see for example Fig. 16 of iMevnet & Maeder] 
12.0031) . The N/He ration thus evolves from the initial value to a 
value corresponding to the CNO equilibrium. Stars which have 
not yet reached this equilibrium will show intermediate values 
of N/He. Focusing on the left part of Fig. [8] corresponding to 
the case of solar metallicity tracks, we see again a clear dif- 
ference between the WNLh and the O stars: the latter appear 
to lie on evolutionary tracks where the N/He ration is still ris- 
ing, while the former are all in the region where N/He is con- 
stant. This confirms our previous finding that the WNLh stars are 
core H burning stars which have reached the CNO equilibrium. 
On the contrary, the O stars are still in the process of reaching 
this equilibrium. An important discrepancy between evolution- 
ary tracks and observed stars is the very large values of N/He in 
the Arches WNLh stars compared to the tracks: for some stars, 
no track seems to be able to reproduce their position. Even tracks 
for stars more massive than 120 M (the most massive star for 
which evolutionary tracks exist) would not help, since all the 
tracks between 60 and 120 M seem to predict about the same 
amount of N/He at equilibrium. More massive stars would show 
the same values. If we turn to the right part of Fig. |8]where the 
twice solar metallicity tracks are shown together with the derived 
stellar properties of the Arches stars, the situation dramatically 
improves. All stars can now be represented by the theoretical 
tracks, even the WN7-9h stars with the largest N/He ratios. The 
ones with the lowest ratios now appear not to have reached com- 
pletely the CNO equilibrium. The reason for the largest theo- 
retical N/He ratio at twice solar metallicity is that between Z 
and Z = 2 x Z the initial He fraction barely changes, while the 
N content is much larger. The initial N/He ratio is thus larger. 
Does that mean that the Arches cluster metallicity is super so- 
lar? Although Fig. [8] makes a good case for it, we have seen 
that the interpretation was a bit different for Fig. |4]for which a 
super solar metallicity did not fully explain the observed trend. 
The question of the Arches cluster metallicity will be debated in 
Sect.0 

A final plot which combines the two previous ones is shown 
in Fig. [6] It displays the abundance ratio N/C as a function of T e ff . 
Here again, this ratio reaches a maximum value when the stars 
is at the CNO equilibrium. The WNLh stars and the O stars are 
clearly separated, the former being more N rich-C poor than the 
latter. Note however the existence of some overlap, two WNLh 
stars having N/C ratios similar to O stars. We note that some 
WN7-9h stars still show N/C ratios larger than predicted. Using 
evolutionary tracks with Z = 2 x Z does not really improve 
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Fig. 5. Logarithm of the ratio of nitrogen to helium mass fraction in the Geneva evolutionary models of Meynet & Mae derl (120051) 
and as derived in the stars analyzed in this work (symbols). Left: solar metallicity evolutionary tracks. Right: twice solar metallicity 
tracks. 



the situation in this case, since the initial N and C abundances 
are changed in a very similar way when Z is increased, so that 
the initial N/C ratio does not vary a lot. Finally, it is interesting 
to note that the scatter of the N/C ratio seen in Fig. [6] among 
WN7-9h stars is real. An observational demonstration is made 
in Fig. [7] where we see that the two stars F6 and F7 have very 
similar spectra (and thus similar parameters, see Table [2]) but 
different Civ 2.070-2.084 fim and Nm 2.247, 2.251 fim lines. 
This means that stars at the same position in the HR diagram 
can be in slightly different states of chemical evolution. 

From this analysis, we can safely conclude that the WN7-9h 
stars in the Arches are core H burning objects showing products 
of the CNO equilibrium at their surface. They are clearly distinct 
from the less evolved O supergiants. 

6.2. Winds 

In Fig. ?? we show the mass loss rates of the Arches stars studied 
here as a function of luminosity (left). It is clear that the WNLh 
stars have stronger winds than the O stars. Indeed, although mass 
loss increases with lumi nosity as predict ed by the theory of ra- 
diation driven winds (Casto r et al.lll975t) . there is a separation 
between the O and WNLh stars for the luminosity range 5.7-6.2 
in which both types of stars are found. This is another indication 
that WN7-9h stars are more evolved than O stars. It is also very 
interesting to note that the extreme supergiants (04-6If + stars) 
seem to have mass loss rates intermediate between normal O su- 
pergiants and WN7-9h stars. As we will see in Sect. 16.31 this 
is an indication of an evolutionary link between early O super- 
giants and WNLh stars. 

On the right part of Fig. ??, we show the clumping corrected 
mass loss rates (M/ ^ff) as a function of luminosity for our sam- 
ple stars as well as for comparison objects. We see that quali- 
tatively, the mass loss rates we derive are consistent with pre- 
vious studies of Galactic stars. Given the scatter in M among 




4.7 4.6 4.5 4.4 4.3 



log (T e „) 

Fig. 6. Logarithm of the ratio of nitroge n to carbon mass fraction 
in the Geneva evolutionary models of Meynet & Maeder (2005) 
and as derived in the stars analyzed in this work (symbols). 

the various types of stars, it is not possible to make a quanti- 
tative comparison. Instead, the so-called modified wind momen- 
tum - luminosity relation, or WLR, is a better tool. The modified 
wind momentum, Mvoo V/? (R being the stellar radius), is ex- 
pected on theoretical gr ounds to depend only on luminosity (e.g. 
iKudritzki & PulslfeOOOh . Fig.|9]shows the relation for the Arches 
cluster stars. On average, all stars follow the qualitative increase 
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Fig. 8. Left: Mass loss rates as a function of luminosity for the Arches cluster stars. Right: clumping corrected mass loss rates 
for t he Arches clu s ter sta r s (filled symbo l s) as well as compa r ison o bjects (open symbols): e arly type (extreme) supergiant s 
from iHerrero et al.1 (120021) : iRepolust et ail d2005l) : iBouret et all (120051) and WN6-9h stars from ICrowther & Bohannanl d 1997b : 
ICrowther & Dessartl (119981) 




modified wind momenta and luminosity. More specifically, we 
have: 



Fig. 7. Comparison between the spectra of stars F6 and F7 show- 
ing that stars with similar main lines and thus similar parameters 
(T e ff, log M) can have different C and N abundances, as re- 
vealed by the different strengths of the C iv 2.070-2.084 ymi and 
Nm 2.247, 2.251 ym lines. 



of the modified wind momentum as a function of luminosity. It 
is very important to note that this is also true to a large extent 
for the WN7-9h stars: they clearly show a correlation between 



logMv ra Vft = 



21.35(±2.54)+ 1.34(±0.44)log^ OI 
23.33(±1.28)+1.08(±0.21)log£ WNLh 



(1) 



The existence of a WLR for WN7-9h stars is a strong in- 
dication that radiative acceleration plays an important in driv- 
ing their winds. This is a fundamental result, since the question 
of whether or not radiative acceleration is efficient enough to 
produce the large WR stars outflows is not entirely settled yet 
(see Crowther at al. 2007 for a recent revie w). Recent theoretical 
simula tions by lGrafener & Ha mann (2005) and lVink & de K oter 
(2005) convincingly indicate that radiative acceleration might 
explain the large mass loss rates of WR stars. Here, we pro- 
vide an observational evidence that at least qualitatively, the O 
and WNLh stars winds in the Arches cluster rely on the same 
physics. Another argument in favor of line driving for the winds 
of the WN7-9h stars comes from the values of 77 = Mx v x /(L/c) 
listed in the last column of Table[2] We see that for these stars, it 
is close to 1, meaning that radiation alone is in principle able to 
drive the wind, even in the single scattering limit. 

The WLR followed by the WN7-9h stars is systematically 
shifted towards higher values (by ~ 0.4 dex) compared to the O 
stars relation. This is not surprising since such an effect is also 
seen among O stars with d ifferent luminosity classes as well as 
betwe en O, B and A stars ( IRepolust et alj|2004t iKudritzki et al.l 
1999). This is usually interpreted as a change in the number 
of lines effectively dri ving the acceleration (see discussion in 
IKudritzki & Puil2000l) . We also note that the WLR we derive is 
flatter for WN7-9h stars: the slope is 1.34 for O stars, and 1.08 
for WN7-9h stars. Given the errors, we cannot exclude however 
that the slopes are similar. Note that if we calculate the slope of 
the WLR regardless of the spectral types, i.e. including all stars 
of our sample, we get a value of 2.00±0.19. 
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In Fig. [9] we have also plotted the theoretical relation of 
IVink et alj (BoOOl) for O stars. Its slope is 1 .83+0.04. We see that 
given our limited sample and the errors, our WLR slope for O 
stars is rather similar. The absolute position of the WLR is how- 
ever different. This is a known effect usually attrib uted to the 
use o f clumping in our models, while the models of IVink et al.l 
(2000) are homogeneous. Note however that hydrodynamical 
confirmations of this effect are still lacking. Including clump- 
ing in atmosphere models leads to systematically lower values 
of M compared to studies with homogeneous winds, the differ- 
ence being Jf with / the clumping factor. In our case, we use 
/ = 0. 1 so that we can expect a shift of 0.5 dex in the WLR of 
O stars (see dotted line in Fig . [9j> . We see that if we correct our 
derived relation by this am ount, we end up slightly above the re- 
lation of IVink et all d2000l) . by a factor ~ 0.2 dex (the difference 
between their theoretical relation and ours is ~ 0.2-0.3 dex de- 
pending on the luminosity). How can we explain this remaining 
discrepancy? One might argue that our value of the clumping 
factor, assumed to be 0.1, is not approp riate for O s tars. O ur re- 
lation would be consistent with that of Vin k et al. I (120001) if we 
had chosen / ~ 0.2 - 0.3, a value that we cannot discard from 
our modeling. As previously recalled, we note h owever that re- 
cent studies indicate smaller / for O supergia nts dCrowther et alj 
l2002at iHillier et al.ll2003t iBouret et al.ll2005l) . Another interest- 
ing possibility to explain the difference between our and Vink et 
al.'s WLR is that we see the effects of high metal con tent. Mass 
loss r ates of O stars are expected to scale as Z 85 dVink et alj 
|2001|). If the difference we observe between the theoretical rela- 
tion o flVink et al.l ([2000) and our derived relation for O stars was 
due to such an effect, this would mean that the metallicity of the 
Arches stars should be Z/Z Q = 10"* ~ 1.7. This is an intrigu- 
ing possibility that we will discuss furt her in Sect. [ 7] Fina lly, we 
should mention that the prediction of IVink et al.l d2000t) might 
not be correct. However, several recent stu dies seem to confirm 
its validity, at least at high luminosities (Markovaetal][2004; 
iMokiem et all2006l) . 

To conclude this section on the wind properties of the Arches 
stars, we summarize in Table [4] the mass loss rates of the stars 
analyzed in the pres ent study which h ave been observed at ra- 
dio wavelengths by lLang et all (120051) . From these detections, 
the authors der ive values of the mass lo ss rates using the stan- 
dard relation of iWright & Barlowl(ll975l) . Since this relation im- 
plies that stellar wind s are homogeneou s, we have to compare 
the mass loss rates of Lang et al. (20051) to our clumping cor- 
rected M (i.e. -^=). lLang et all ((2005) also use a different dis- 

tance (8.0kpc instead of 7.62kpc) and assume p - 2.0 to derive 
their mass loss rates. We also hav e to correct for thi s. The ini- 
tial and corrected values of M of lLang et ail ([2005) are given 
in columns 3 and 4 of Table [4] For stars 3, 5, 8 and 18, we see 
that both our and the radio determination are consistent within 
a factor of 2. For the remaining stars (1, 2, 4, 6), the differ- 
ences can be as large as almost an order of magnitude. Star 
F6 is suspected to be the c ounterpart of the X-ray source A IN 
dLaw & Yu sef-Zadeh 2004), indicating tha t it might be a n on- 
thermal emitter (colliding winds). Besides. lLang et all d2005h re- 
port that star F6 is variable. Hence, differences between our de- 
termination and the radio mass loss rate is not surprising. For the 
three remaining stars (Fl, F2 and F4), we have no explanation 
so far of the disagreement. One might invoke binarity (although 
our spectra exclude the presence of a spectroscopic companion) 
or crowding which may affect the radio determination (the reso- 
lution - 0.42" x 0.17" - being lower than in our SINFONI data 
- 0.20"). 
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Fig. 9. Modified wind momentum - luminosity for the Arches 
cluster stars. The solid line is the theoretical prediction from 
IVink et all d2000t) . The red dashed line is the linear fit to the O 
stars, excluding the 04-6If + supergiants. The blue dot dashed 
line is the fit to the WN7-9h stars. The dotted line is our WLR 
for O stars if we use unclumped mass loss rates. See text for 
discussion. 



Excluding star F6 (suspected binary), it is interesting to note 
that the radio mass loss rates are systematically lower than the 
IR unclumped on es (exception: F18). This effect was noted by 
Fige r et aTTd2002l) for one star (F8). In our case, he difference is 
0.36 dex, with a rather large dispersion (0.32). Since the ratio of 
unclumped mass loss rates is virtually similar to the ratio of wind 
densities (recall that the M diagnostics are in fact density indica- 
tors), this difference might tell us something about the different 
clumping factors in the IR and radio emitting regions. Density 
is oc Ml y/, where M is the "true" mass loss rate. Hence, the 
ratio of radio to IR unclumped mass loss rates is oc *>Jf IR / f mc i io 
where //« (f ra dio) is the clumping factor in the IR (radio) emit- 
ting region. With the observed trend, it seems that clumping is 
stronger in the IR than in the radio emitting region (by a factor 
of ~ 5). This is in agreem ent (qualita t ively and quantitatively) 
with the recent findings of Puis et al. (2006J ) and with the the- 
oretical predictions of iRunacres & Owockil d2002l) . This might 
partly explain the systematic differences seen in Table [4] 



6.3. Nature and evolution of the most luminous stars 

The discussion in the last sections has lead to a clear picture 
for the nature of the WN7-9h stars in the Arches cluster: they 
are very massive (60 < M < 120M Q J post main sequence ob- 
jects still in the H burning phase and have reached the CNO 
equilibrium. They clearly separate from the rest of the O stars 
studied here. As such, they are very reminiscent of the H rich 
WN stars in the core of the NGC 3603 and R136 clusters 
dCrowther & Dessart|[T998T) . 
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Table 4. Clumping corrected mass loss rates for the st ars in com- 
mon b etween the present sample and the sample of lLang et al.l 
(120051) . 



Star log -y= log M (Lang et al.) log M (Lang et al.) corrected 



1 -4.20 -4.72 -4.80 

2 -4.22 -4.72 -4.66 

3 -4.10 -4.37 -4.47 

4 -3.85 -4.72 -4.64 

5 -4.14 -4.43 -4.44 

6 -4.12 -3.65 -3.66 
8 -4.00 -4.34 -4.27 
18 -4.85 -4.72 -4.62 



From the anal ysis of several Galactic WN stars, 
ICrowther et alj d 19951) built the following evolutionary se- 
quence for stars more massive than 60 M : 

O Of WNL+abs WN7 WNE) WC -> SN 

where WNL ( WNE) stands for W N late (early), i.e. WN6-9 
(WN2-5) stars. iLanger et alj dl994l) describe a similar sequence 
except that 1) they include a Luminous Blue Variable phase af- 
ter the H-rich (=WNL+abs) phase, and 2) they do not explicitely 
list WN7 types in their scenario but replace them by the term H- 
poor WN. However, they tentatively identify H-rich WN stars as 
core H burning objects. They also note that these objects are the 
most luminous WN stars. 

The Arches WNLh stars nicely fit the global picture drawn 
in these scenarios: they are very luminous, young, H rich objects 
clearly still in the H burning phase and with masses in excess of 
60 M Q . Besides, our detailed analysis of the C and N abundance 
patterns of these stars quantitatively strengthens the conclusion 
that they are in a relatively early evolutionary state. The direct 
link between O and WNLh stars is also confirmed by our analy- 
sis. The O supergiants of our sample all appear to be less evolved 
than the WNLh stars (see Fig. [3j 13 [8] and©. And the two 04- 
6If + extreme supergiants have properties intermediate between 
normal supergiants and WNLh stars: this is best seen in Fig. [8] 
where they bridge the two latter classes of objects. Note how- 
ever that we cannot state that there is a direct link between spe- 
cific sub-classes of stars in the Arches cluster, namely between 
04-6 supergiants and WNLh stars: Fig. [2] reveals that the for- 
mer probably have lower initial masses than the latter. But we 
can safely conclude that WNLh stars have early O stars as pro- 
genitors, since they should evolve from more massive (and con- 
sequently hotter) O stars than the O supergiants of our sample. 
In practice, the Arches W NLh stars should be the desce ndant of 
02-4 stars. We note that ICrowther & Bohannanl d 19971) argued 
for a direct link between 08If and WN9ha stars. However, this 
was based on the study of only 3 stars, compared to 28 in the 
present paper. The WN9ha of their study had in addition a lower 
luminosity and temperature than the Arches WN7-9h stars, and 
might thus be a different, initially less massive type of WN9h 
star. Conseq uently, we do not think that our results are in contra- 
diction with Crowt her & Bohannanl dl997l) . but rather that they 
refer to different kinds of stars. 

In conclusion, our findings strongly support the scenario ac- 
cording to which, in the Arches cluster, the most massive O stars 
evolve into extreme supergiants and then into H-rich WNL stars. 



7. Stellar metallicity 

In the previous sections we have seen that several elements 
pointed towards a super solar metallicity for the Arches cluster 
stars: the N enrichment might be too large to be accounted for 
by solar metallicity evolutionary models, and the winds might be 
stronger than expected for a solar compositi on. This is somewha t 
in contradiction with the recent results of Naiar ro et al.l (2004) 
who favored a solar metallicity for the five stars they analyzed. 
Their determination was based on the interesting finding that the 
N mass fraction reaches a maximum in evolutionary tracks when 
the star is in the WN phase. This maximum does not depend on 
the initial mass, but is sensitive to the initial global metallicity. 
In practice, comparing the mass fraction of a sampl e of WN stars 
to such tracks should then constrain the metallicity. Naiarr o~et al.l 
(2004) used the three WN stars they analyzed to make such an 
estimate. They found that a solar metallicity was preferred. In 
Table. [3] we have shown that for these three stars, we find similar 
X(N) for one (F8), and slightly larger values than Naiarro et al. 
(2004) for the other two (F3 and F4). Fig.[10]is the figure Najarro 
et al. used to estimate Z in the Arches cluster, but now using 
the 13 WN stars of our sample: the shaded area corresponds to 
the range of X(N) covered by these 13 stars. From this, we see 
that a wide range of metallicity is possible. If real, this can be 
attributed to two factors. First, there may be a scatter in the ini- 
tial metallicity of the Arches stars. However, a difference of a 
factor of 2 seems quite large. The second possibility is that the 
WN stars might not all have reached the phase of their evolution 
were X(N) is maximum. In that sense, Fig. [TOl provides only a 
lower limit on the metallicity. Note that this effect most likely 
influences the results when a small number of stars is used. One 
might also wonder whether the scatter we see is not purely sta- 
tistical. If we assume it is the case, we can run a x 2 analysis to 
find the preferred metallicity, using 

X 2 (Z) = V ( *'-*7 (Z))2 ( 2 ) 

where Xj are the individual N mass fraction, <x, the associated un- 
certainties, and X max {Z) the maximum N mass fraction reached 
at a given metallicity. The evolution of x 1 In with Z is shown in 
Fig. [TT](« being the number of free parameters, equal to 12 in 
our case - 13 stars included and one fitted parameter). We see 
that there is a clear minimum for Z=l. 3-1.4 Z , and that a twice 
solar metallicity is clearly excluded. The minimum x 1 l n is 1-88, 
so a unique metallicity is not completely satisfactory to explain 
the distribution of X{N) max . One has to keep in mind that this 
analysis is also correct only if the evolutionary models correctly 
predict the evolution of surface abundances. 

Altoget her, our resul t s ten d to favor a slightly super solar 
metallicity. iNaiarro et al.l (120041) concluded that Z=Z was pre- 
ferred. The difference is likely due to the use of different evolu- 
tionary tracks: we use the recent tracks from lMevnet & Maederl 
(200 5j), while Najar ro et al. use the obsolete tracks from 
Schal ler et al.l (1992). Since for a given metallicity, these old 
tracks have a larger maximum X(N), a lower metallicity is de- 
rived for a given range of observed X(N) In our sample, the 
average value of X(N) in WNLh stars is 1 .70(±0.60)%. The three 
stars of Najarro et al. have on average X(N)=1.57(+0.15)%. 
Hence, within the errors, the values are similar. The different de- 
rived metallicity we find is thus likely due to the use of different 
evolutionary tracks. 

3 Note that the maximum value of X(N) is independent of the rotation 
rate of the star. 
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An important comment to make is that the approach de- 
scribed above is only valid if the CNO abundances and the global 
metallicity Z scale similarly, or stated differently if all the met- 
als have the same relative overabundance compared to the so- 
lar composition. Let us assume we have a model in which the 
initial CNO content is larger than solar, while all the other met- 
als have solar abundances. The global metallicity will then be 
slightly above solar. Let us consider a second model in which 
the global metallicity is the same as in the first one, but now all 
the elements have the same abundance excess relative to the so- 
lar composition. In the two models, the global metallicity is the 
same, but the first one will produce a much larger N mass frac- 
tion. If we use this mass fraction to asses the global metallicity 
using the method presented above, then we will overestimate the 
global metallicity of the star. 

Studies of abundances of individual elements not affected by 
stellar evolution would help to better constrain the metallicity of 
the Arches cluster stars. Unfortunately, very few lines are avail- 
able for such a purpose. The only one which could be of help 
is the Si iv line at 2.428 yum. Interestingly, it turns out that to 
correctly reproduce it in all the stars we have studied so far, we 
need a silicon abundance between solar and twice solar. Given 
that the line is located at the end of the K band where the S/N 
degrades rapidly and where the spectrum normalization is less 
straightforward than at shorter wavelengths, we claim that this 
is only an indication that individual abundances of light metals 
might be super-solar. It is nonetheless interesting that this would 
be consistent with the various indications gathered so far. If we 
assume that the Fe content is abou t solar as studies of red super- 
giants in the central cluster show dCarr et al.l2000tlRamfrez et alj 
2000), it might well mean that there is an excess of light metals 
relative to heavier metals in the Arches cluster. Such a conclu- 
sion would thus weak en the results of the Z determination by the 
method presented by Naiarr oet alj d2004l) . But this would be a 
very attractive possibility, since this could indicate a larger a/Fe 
abundance ratio, which in turn is usually interpreted as the im- 
print of a top-heavy IMF. A recent study of stellar abundances 
in cool luminous stars bv lCunha et alj d2007f) concluded that O 
and Ca (two a elements) were overabundant compared to Fe in 
the central cluster. Our suggestion of a super-solar Si abundance 
in the Arches cluster is consistent with their findings. 

In conclusion, we tentatively suggest that in the Arches clus- 
ter, the lightest elements most likely have a super-solar abun- 
dance while the iron peak elements have a solar metallicity. 

8. Conclusions 

We have presented a quantitative study of the most massive stars 
in the Arches cluster. K-band spectra have been obtained with 
SINFONI on the VLT. A detailed spectral classification has re- 
vealed the presence of WN7-9h stars as well as O supergiants, 
including two extreme OIf + stars. We have quantified the main 
stellar and wind parameters of 28 stars using the atmosphere 
code CMFGEN. The main results of our study are: 



0.03 



0.02 



0.01 




Age [Myr] 

Fig. 10. N mass fraction as a function of age in evolution- 
ary models for solar (solid line) and twice solar (dotted line) 
m etallicity. The Geneva evo lutionary tracks including rotation 
of iMeynet & Maederl d2005l) are used. Tracks for M= 120, 85, 60 
and 40 M are shown from from left to right. The shaded area in- 
dicates the range of X(N) covered by the WN7-9h stars analyzed 
in the present study. 




- the massive star population of the Arches cluster is 2—4 Myr 
old. Although marginal, there seems to be the trend that the 
most massive stars are also the youngest: WNLh stars are 
2-3 Myr old, while less massive O supergiants might have 
an age of up to 4 Myr. Initial masses as large as 120 M are 
derived for the WNLh stars from the HR diagram. 

- The WN7-9h stars are identified as core H burning stars 
which show chemical enrichment typical of the CNO equi- 
librium: they still contain a significant amount of hydrogen 



Fig. 11. Determination of metallicity from the N mass fraction of 
WNLh stars, x 1 is defined in Eq.|2] x 2 l n 1S shown as a function of 
metallicity. n is the number of free parameters (12 is the present 
case). If the scatter in X(N) is purely statistical, then a slightly 
super-solar metallicity (Z=1.3~1.4Z©) is preferred. 
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and show both N enhancement and C depletion. They have 
supergiants of spectral types earlier than 04-6 as progeni- 
tors. 

- The level of N enrichment suggests either a super-solar ini- 
tial content and/or a too low efficiency of N enrichment in the 
evolutionary models. The indication that Si might be over- 
abundant by a factor 2 compared to the solar abundance ar- 
gues in favor of a super-solar metallicity at least for the light- 
est metals. 

- The properties of the Arches massive stars argue in favor of 
the evolutionary scenario of ICrowther et al.l (I 19951) for the 
most massive stars: O — > Of — > WNL+abs — > WN7 

- The winds of the WN7-9h stars follow a well defined modi- 
fied wind momentum - luminosity relation. This is a strong 
indication that they are radiatively driven. It also seems that 
the winds are less clumped in the radio continuum emitting 
regio n than in the nea r-IR line emitting region, in agreement 
with lPuls etail d2006l) . 

In order to test the indication that the most massive stars are 
the last to form, intermediate mass stars must be analyzed. The 
question of the metallicity in the Galactic Center, and the Arches 
cluster in particular, is clearly not answered yet. New observa- 
tions/analysis of cooler stars showing a larger number of metallic 
lines are needed. We will present such a study in a forthcoming 
paper. 
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Appendix A: Best fits 

In this Section, we gather the figures showing the comparison 
of our best fit models with the observed spectra of our program 
stars. 
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Fig. A.l. Best fits (red dot-dashed lines) of the observed K-band spectra (black solid line). 
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Fig. A.2. Best fits (red dot-dashed lines) of the observed K-band spectra (black solid line). 
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Fig. A.3. Best fits (red dot-dashed lines) of the observed K-band spectra (black solid line). 
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